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Abstract 

Recently a one-dimensional closed ladder of Josephson junctions has been 
studied [1] within a twisted conformal field theory (CFT) approach [2, 3] and 
shown to develop the phenomenon of fiux fractionalization [4]. That led us to 
predict the emergence of a topological order in such a system [5]. In this letter 
we analyze the ground states and the topological properties of fully frustrated 
Josephson junction arrays (JJA) arranged in a Corbino disk geometry for a 
variety of boundary conditions. In particular minimal configurations of fully 
frustrated JJA are considered and shown to exhibit the properties needed in 
order to build up a solid state qubit, protected from decoherence. The stabihty 
and transformation properties of the ground states of the JJA under adiabatic 
magnetic fiux changes are analyzed in detail in order to provide a tool for the 
manipulation of the proposed qubit. 
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1 Introduction 



Today the physical reahzation of a quantum computer represents a very hard task because 
of the hmits imposed by dccohcrence and because interactions between qubits cannot be 
well controlled. A lot of solid state qubit implementations have been proposed in the last 
ten years; among them, systems based on Josephson junctions [6] are promising because 
of the well assessed fabrication technology, which could allow a high degree of scalabil- 
ity. But anyway the path towards the physical realization of a fully integrated quantum 
computer still remains a remarkable challenge. Recently topologically ordered quantum 
systems have been proposed as physical analogues of quantum error- correcting codes [7] 
starting from Kitaev seminal work [8]: all such systems share a built in protection from 
decoherence. The idea of Kitaev involves a protected subspace created by a topological 
degeneracy of the ground state: such a degeneracy is typically due to a conservation law 
such as the conservation of the parity of the number of particles along some long contour. 
The underlying idea in Kitaev work is the notion of topological order, that is a new kind 
of order which cannot be due to a spontaneous symmetry breakdown. In such a case 
the system lacks of local order parameters but displays a weak form of order, which is 
sensitive to the topology of the underlying two dimensional manifold. Such a concept 
was first introduced in order to describe the ground state of a quantum Hall fluid [9] but 
today it is of much more general interest [10] [11]. Two features of topological order are 
very striking: fractionally charged quasiparticles and a ground state degeneracy depend- 
ing on the topology of the underlying manifold, which is lifted by quasiparticles tunneling 
processes. In general a system is in a topological phase if its low-energy, long-distance 
effective field theory is a topological quantum field theory that is, if all of its physical 
correlation functions are topologically invariant up to corrections of the form e~rat tem- 
perature T for some nonzero energy gap A. More recently superconductors have been 
proposed [12] in which superconductivity arises from a topological mechanism rather than 
from a Ginzburg-Landau paradigm: the key feature is a mapping on an effective Chern- 
Simons gauge theory, which turns out to be exact in the case of JJA and frustrated 
JJA [13]. Non-Abelian quantum Hall states, in particular the t-' = \ one, appear also 
very promising in that the quasiparticle excitations are non-Abelian anyons [14] obeying 
non-Abelian braiding statistics: in such a case quantum information is codified in states 
with multiple quasiparticles, which show a topological degeneracy. Furthermore topo- 
logical phases have been recently recognized in chiral p-wave superconductors [15], such 
as Sr2Ru04, and ultra-cold atomic gases [16]. Large and small size Josephson junction 
arrays of special geometry [17] [18] have been proposed as well, which share the property 
that, in the classical limit for the local superconducting variables, the ground state is 
highly degenerate. The residual quantum processes within such a low energy subspace 
lift the classical degeneracy in favor of macroscopic coherent superpositions of classical 
ground states [18]. The protected degeneracy in all such systems emerges as a natural 
property of the lattice Chern-Simons gauge theories which describe them [18]. In general, 
if a physical system has topological degrees of freedom that are insensitive to local pertur- 
bations (that is noise), then information contained in those degrees of freedom would be 
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automatically protected against errors caused by local interactions with the environment 
[8]. The procedure implemented in all such realizations of protected JJA based qubits 
runs as follows. A quantum system is required with 2^ quantum states [K being the num- 
ber of big openings in the Josephson systems under study) which are degenerate in the 
absence of external perturbations and are robust against local random fluctuations. This 
means that the Hilbert space should contain a 2^-dimensional subspace characterized by 
the crucial property that any local operator O has only state-independent diagonal matrix 
elements up to vanishingly small corrections: {n\0\m) = OoSmn + o [exp {—L)], L being 
the system size. So a possible answer to such a highly non trivial requirement could be 
a system with a protected subspace built up by a topological degeneracy of the ground 
state [8] . An alternative procedure could be to exhibit a low-energy effective field theory 
for the system under study which is a topological one and whose vacua are topologically 
degenerate and robust against noise. Following such an approach, we were able to predict 
the emergence of topological order in a fully frustrated Josephson junction ladder (JJL) 
with Mobius boundary conditions [5, 4] and to build up a protected subspace with 2^ 
quantum states, K — 1, which could be identified with the two states of a "protected" 
qubit [1] . The low-energy effective field theory which we set up is a twisted conformal field 
theory, the Twisted Model (TM) [2, 3]: it accounts very well for the topological proper- 
ties of the system under study [5, 4]. This finding is crucial and relies on the well known 
relation between Chern-Simons gauge theories in (2 -|- l)-dimensions and 2-dimensional 
CFT [19]. 

The aim of this letter is to analyze the ground states and the topological properties of 
fully frustrated Josephson junction arrays (JJA) arranged in a Corbino disk geometry for 
a variety of boundary conditions, employing a twisted CFT approach [2, 3] which has been 
successfully applied to quantum Hall systems in the presence of impurities or defects [20, 
21, 22] and to the study of the phase diagram of the fully frustrated XY model {FFXY) 
on a square lattice [23]. In this way we extend to a more general and more involved 
system the results obtained in the simplest JJL case [5, 4, 1] and exhibit the minimal 
configuration and properties needed in order to build up a solid state qubit, protected 
from decoherence. We analyze in detail the stability and transformation properties of the 
ground state wave functions under adiabatic magnetic fiux changes in order to provide a 
tool for the manipulation of the proposed qubit [24] . 

The letter is organized as follows. 

In Section 2 we briefiy discuss all the relevant phenomenology of 2-dimensional fully 
frustrated JJA with an emphasis on the underlying physical model, the FFXY model 
on a square lattice [25]. The analysis will focus on some models which share the same 
t/ (1) (g) Z2 degenerate ground state [26] and are believed to be in the same universality 
class. 

In Section 3 we recall some aspects of the m-reduction procedure, in particular we 
show how the m = 2, p = case gives rise to the required U (1) ^ Z2 mixed symmetry 
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[23] of the FFXY model. Then we focus on the discrete version of such a procedure. In 
such a framework we give the whole primary field content of the theory on the plane. 

In Section 4 we properly define a JJA on the torus topology, which is equivalent to the 
Corbino disk geometry with coinciding boundary, and introduce the magnetic translation 

operators. We build up the ground state wave functions in order to provide a physical 
identification of the characters of our CFT, the TM, as the components of the "center of 
charge" for such wave functions. Then we describe the general topological structure of 
the vacua on the torus corresponding to all the possible boundary conditions. 

In Section 5 we study the stability and transformation properties of the four ground 
states of the JJA, arranged in the closed geometry under adiabatic magnetic flux changes 
through the central hole of the Corbino disk. That allows us to identify the two states 
of a possible flux qubit [27], protected from decoherence, and to provide a tool for its 
manipulation. 

In Section 6 some comments and outlooks are given while in the Appendix the explicit 
expressions of the TM characters on the torus are presented in detail. 

2 Fully frustrated Josephson junction arrays 

Josephson junction arrays (JJA) are prototype systems to investigate a variety of phase 
transitions induced by thermal or quantum fluctuations [28]. If the superconducting 
islands are of submicron size, quantum fluctuations play a crucial role and drive the zero 
temperature Superconductor-Insulator (SI) phase transition [29]. There are two energy 
scales in JJA: the Josephson energy J, which is associated to the tunneling of Cooper 
pairs between neighboring islands, and the charging energy U , which is the energy needed 
in order to add an extra Cooper pair on a neutral island. A flnite U is responsible 
of quantum fluctuations of the phases {</?} of the superconducting order parameter on 
each island. In the regime J ^ U, named the classical case, the fluctuations of the 
phases are small, the system is globally coherent and superconducting. Conversely, in 
the opposite regime J -C t/, strong quantum phase fluctuations prevent the array from 
reaching long-range phase coherence and make it a Mott insulator. Magnetic frustration 
can be introduced in a JJA by applying a magnetic fleld transversal to the plane of the 
array [30] [31]; it can be deflned as / = where $ is the magnetic flux threading each 
plaquette and $o = |f is the superconducting flux quantum. The external magnetic fleld 
induces vortices in the array and, in the case of rational frustration / = |, the ground 
state shows a checkerboard configuration of vortices on a. q x q elementary supercell. In 
the case of full frustration, i.e. / = |, of interest to us here, there are two degenerate 
ground states built of a vortex lattice with a 2 x 2 elementary supercell; the corresponding 
current flows either clockwise or anticlockwise in each plaquette, giving rise to a chiral 
phase conflguration. Quantum fluctuations affect the superconducting regime but do not 



4 



destroy the checkerboard structure of the ground state [32]. All the relevant physics of 
the JJA is captured by the quantum phase model (QPM), whose Hamiltonian is: 

Hqpm = XI (^i - n^) Uij {rij - n^;) - J ^ cos {(pi - (pj - Aj) , (1) 

(ij) 

where and (pi are canonically conjugated variables defined on the sites and satisfying the 
commutation relations rij] = 2ei6ij, the second sum is over nearest neighbors, J > is 
the Josephson energy, the matrix Uij = Ae^C'^^ describes the Coulomb interaction (Qj is 
the capacitance matrix), the external voltage enters through the induced charge eux 
and fixes the average charge on each island and Aij — A-dl is the line integral along 
the bond between adjacent sites i and j. We consider the case where the bond variables 
Aij are fixed, uniformly quenched, out of equilibrium with the site variables and satisfy 
the condition = 27r/; here the sum is over each set of bonds of an elementary 

plaquette and / is the strength of frustration. In the following we will focus on the limit 
J ^ U, where the JJA in the presence of an external magnetic field transversal to the 
lattice plane is a physical realization of the fully frustrated XY model on the square 
lattice. Let us then focus on such a model, which is described by the action: 

if = - J X COS {tpi - ipj - Aij) . (2) 

(iJ) 

We assume that the local magnetic field in Eq. (2) is equal to the uniform applied field; 
such an approximation is more valid the smaller is the sample size L compared with the 
transverse penetration depth Aj^. In the case under study, / = |, such a model has a 
continuous U{1) symmetry associated with the rotation of spins and an extra discrete 
Z2 symmetry, as it has been shown analyzing the degeneracy of the ground state [25]. 
Choosing the Landau gauge, such that the vector potential vanishes on all horizontal 
bonds and on alternating vertical bonds, we get a lattice where each plaquette displays one 
antiferromagnetic and three ferromagnetic bonds. Such a choice corresponds to switching 
the sign of the interaction term in Eq. (2) and is closely related to the presence of two 
ground states with opposite chiralities, the first one invariant under shifts by two lattice 
spacings and the second one invariant under shifts by one lattice spacing. The action (2) 
can be cast into a form where both the U{1) and Z2 symmetries are manifest, through 
the Villain approximation [33]. In this way the spin- wave and the vortex contributions 
can be separated. Furthermore, by integrating out the spin waves, the resulting vortex 
contribution can be rewritten as a fractionally charged Coulomb gas (CG) defined on the 
dual lattice [34]: 

H^-jJ2(m(r) + f)G (r, r') (m (r') + /) . (3) 

r,r' 

Here we have lim|r-r'|^oo G (r, r') = log |r — r'|+|7r and the neutrality condition J2r (^) + /) — 
X^j. n{r) — must be satisfied. It is now evident that the ground state for / = | consists 
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of an alternating lattice of logarithmically interacting ±| charges and is doubly degener- 
ate. Such a model exhibits two possible phase transitions, an Ising and a vortex-unbinding 
one [35] , and their relative order has been deeply studied in the literature [25] . 

The FFXY model has been studied analyzing other models which have the same 

U{1) (S> Z2 degenerate ground state and are believed to be in the same universality class. 
In particular it can be reformulated in terms of a system of two coupled XY models with 
a symmetry breaking term [36]: 



J2 5]cos((^«(r)-(^»(r')) 
i=l,2 ^^y^ 



+ hY,cos2{(p^'\r)-if^^\r)) . (4) 



The limit /i — >• corresponds to a full decoupling of the fields Lp^'^\ i = 1,2, so giving rise 
to a CFT with central charge c — 2 which describes two independent classical XY models 
(in the continuum limit). In the /i — > 00 limit the two phases ip^^\ i — 1,2 are locked 
[36], i. e. (f^^\r) — (p^'^\r) = nj, j = 1,2; as a consequence the model gains a symmetry 
[/(I) (g) Z2 and its Hamiltonian renormalizes towards a model described by: 

H^H {h^ 00) ^aJ2 {l + SrS^')cos(^ip^^\r) -ip^^\r')^ (5) 

(r,/) 

where (p^^\r) and Sr = costtj = ±1 are planar and Ising spins respectively. In this way 
a model is obtained which is consistent with the required symmetry, the XF-Ising one, 
and whose Hamiltonian has the general form [37] : 

HxY-i = X! (1 + ^'■V) cos (ip^'\r) - ^^'Hr')'^ + Cs,v] ■ (6) 

{r,r') 

In the presence of such a symmetry we can use the m-reduction technique [3] which has 
been successfully applied to a quantum Hall fluid in [2, 3, 20, 21] and to a fully frustrated 
Josephson junction ladder with Mobius boundary conditions in [5, 4, 1]. 

Let us now have a look at the full spectrum of excitations of the FFXY model: 
vortices, domain walls, kinks and antikinks. Vortices are point-like defects such that 
the phase rotates by ±27r in going around them [35]. A domain wall is a topological 
excitation of the double degenerate ground state and it can be defined as a line of links, 
each one separating two plaquettes with the same chirality. So through a domain wall the 
alternating structure (the checkerboard pattern) of the ground state is lost. Kinks and 
antikinks are excitations which live on the domain walls and are described by fractional 
vortices with +1/2 and —1/2 topological charge. As we will see, domain wall excitations 
can be generated from the ground state by closing the lattice. This is our strategy: the 
JJA under study can be closed and arranged in a Corbino disk geometry, which is the 
relevant geometry for the implementation of a protected qubit. Imposing the coincidence 
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between the internal and the external edge we obtain the discretized analogue of a torus. 
In this way we generate a variety of boundary conditions which we can classify by noticing 
that, for the ground state, topologically inequivalent circumstances arise for even or odd 
number of plaquettes along the cycles of the torus. In the even case the end plaquettes 
on the opposite sides of the lattice have opposite chirality, while in the odd case they 
have the same chirality. So the ground state on a lattice maps into the ground state on 
the torus only if the torus has an even number of plaquettes along the two cycles. On 
the other hand a straight domain wall is generated along any cycle of the torus which 
corresponds to opposite sides of the lattice separated by an odd number of plaquettes; all 
the possible cases are illustrated in Figs. 1,...,4 in the so-called minimal configurations. 
Such a behavior has to be taken into account by opportune boundary conditions on the 
fields </7(*\ i = 1, 2 at the edges of the finite lattice. These non trivial boundary conditions 
naturally arise when we implement the m-reduction procedure in the discrete case. So 
the closed geometry gives rise to non-trivial topological properties, which appear deeply 
related to the twofold degeneracy of the ground state and are the source of protection 
from external perturbations. 

The next step is to show how our m-reduction procedure will give rise in a natural 
way to such non trivial topological properties for the closed JJA. The first step will be to 
implement the m-reduction in the discrete case and then to perform the continuum limit 
first on the plane and then on the torus. In this way, by Tising the powerful techniques of 
the CFT, we will be able to propose a qubit device protected from decoherence. 

3 The m-reduction procedure for the JJA: a sum- 
mary 

Here we briefly summarize the main results of our theory, the TM, for the fully frustrated 
JJA [23] . We first construct the bosonic theory on the plane and show that it gives rise 
to a system with the required C/ (1) (g) Z2 symmetry. That allows us to describe the JJA 
excitations in terms of the primary fields Va (z). Then we implement a discrete version of 
the m-reduction procedure in order to clarify the central role played by the closed geometry 
and by the non-trivial boundary conditions in the description of the excitations spectrum. 
The topological structure of the vacua on the torus will be outlined in Section 4 together 
with the corresponding physical configurations of the JJA, the detailed expressions of the 
conformal blocks being reported in the Appendix. 

Let us focus on the m-reduction procedure [3] for the special m — 2 case (see Ref. [2] 
for the general case), since we are interested in a system with U (1) (E) Z2 symmetry and 
choose the "bosonic" theory [23], which well adapts to the description of a system with 
Cooper pairs of electric charge 2e in the presence of non-trivial boundary conditions [20], 
i.e. a fully frustrated JJA. As a result of the 2-reduction procedure [2] [3] we get a c = 2 
orbifold CFT, the TM, whose fields have well defined transformation properties under the 
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discrete (twist) group, which is a symmetry of the TM. Its primary fields content can 
be expressed in terms of a Z2-invariant scalar field X{z)^ given by 

X{z)^\{g^'\z) + Q^^\z)), (7) 

describing the continuous phase sector of the theory, and a twisted field 

<^{z) ^\[Q^'\z) - Q^^\z)) , (8) 

which satisfies the twisted boundary conditions 4>{e^'^z) = —4>{z) [2]. The whole TM 
theory decomposes into a tensor product of two CFTs, a twisted invariant one with c = | 
and the remaining c = | one realized by a Major ana fermion in the twisted sector. In the 
c = I sub-theory the primary fields are composite vertex operators V (z) — {z) ip (z) 
or Vqh (z) — Ux (z) a {z), where 

[/^' (^) = ^ : e^"'^(^) : (9) 

is the vertex of the continuous sector with ai = |, / = 1,...,4 for the SU{2) Cooper 
pairing symmetry used here. Regarding the other component, the highest weight state in 
the isospin sector, it can be classified by the two chiral operators: 

2^/z V / 2^/z V / j^^Q-j 

which correspond to two c = ^ Majorana fermions with Ramond (invariant under the Z2 
twist) or Neveu-Schwartz (Z2 twisted) boundary conditions [2] [3] in a fermionized version 
of the theory. The Ramond fields are the degrees of freedom which survive after the 
tunnelling and the parity symmetry, which exchanges the two Ising fermions, is broken. 
Besides the fields appearing in eq. (10), there are the cf [z) fields, also called the twist 
fields, which appear in the quasi-hole primary fields Vqh {z) ■ The twist fields have non 
local properties and decide also for the non trivial properties of the vacuum state, which 
in fact can be twisted or not in our formalism. Indeed the whole TM theory decomposes 
into a tensor product of two CFTs, a twisted invariant one with c = 3/2 (the Moore- 
Rcad (MR) theory with symmetry U{1) ® Z2) and the remaining c = 1/2 one realized 
by a Majorana fermion in the twisted sector. Such a factorization can be unambiguously 
pointed out on the torus topology as we will show in the following [3]. But now let us 
focus on a discrete version of the procedure just outlined in order to clarify the role of 
the closed geometry and of the non-trivial boundary conditions in the description of the 
full spectrum of excitations of the model [23]. To such an extent let (— L/2,0), (L/2,0), 
(L/2, L), (— L/2, L) be the corners of the square lattice C and assume that the fields (^*^*\ 
i = 1, 2 satisfy the following boundary conditions: 

(^(i)(r) =(^(2)^^^ forre£nx, (11) 
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where x is the x axis. The above boundary conditions allow us to consider the two fields 
(p^^^ and (p^'^^ on the square lattice C as the folding of a single field Q, defined on the 
lattice Co with corners (— L/2, — L), (L/2, — L), (L/2,L), (— L/2,L). More precisely we 
define the field Q as: 

^ ' ~ I for r e (-£) n £o. ^ ' 

We can implement now a discrete version of the m-reduction procedure (m = 2) by 
defining the fields: 

X{r) = ^{Q{r) + Q{-r)), (13) 

$(r) = i(Q(r)-Q(-r)), (14) 

where r G £o- The resemblance with the continuum version of the procedure is evident 
and the fields X and $ are symmetric and antisymmetric with respect to the action of 
the generator g : r ^ —r oi the discrete group Z2. The Hamiltonian in Eq. (4) can be 
rewritten in terms of these fields and, for /i = 0, it becomes: 

H^2A J2 cos (x{r) - A'(r')) cos (|$(r) - $(r')) , (15) 

{ry)ec 

which in the continuum limit corresponds to the action of our TM model: 



A 



d'^x. (16) 



It is worth pointing out that the fields X and $ are scalar fields and so the chiral fields 
defined by Eqs. (7), (8) can be seen as their chiral components. Moreover the group Z2 
is a discrete symmetry group, indeed both H and A are invariant under its action. 



4 JJA on the torus topology: magnetic translations 
and wave functions 

The non perturbative ground state wave functions of the JJA system in the torus topology 
represent coherent states of Cooper pairs on the torus. It can be inferred that such ground 
states can be expressed as correlation functions of the primary fields describing the elemen- 
tary particles (Cooper pairs) of JJA. In the following we will explicitly show that the char- 
acters of the theory are in one to one correspondence with such ground states. In particular 
they describe the components of the "center of charge" for the corresponding ground state 
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wave functions [38] , encoding all the topological properties of the J JA system. To such an 
extent let us define for a single Cooper pair on a torus a x 6 an effective mean-field Hamilto- 

nian of the kind H (x,y) — HQ(x,y)+V {x,y), where HQ(x,y) — —ih^ — 2e~^/c /2m 

is the Hamiltonian in the presence of an uniform magnetic field and V {x, y) is a mean- 
field scalar potential such that V {x, y) — V {x + a,y) = V {x,y + b). It is now possible 
to define the magnetic translations operators S = e^^^"-l^ and T = e^^v^l^ along the two 
cycles A (i. e. the real axis x) and B (i. e. the imaginary axis y) of the torus respectively, 
where: 



= - f = -Z/ia,, Qy = T:y^'-^BX=-%Wy^'-^BX ( 1 7) 

and the gauge choice ~^ (x, y) = {—By, 0) has been made. They satisfy the relations: 

S, n {x, y)] = \f, n {x, y)] =0, cSr = e2'^^*'"'/*or5 , (18) 

where is the magnetic flux threading the torus surface, and their action on the wave 
functions can be defined as: 



S(p {x, y) = (p{x + a,y) , Tip {x, y) = (P'-^^Bbx/'i>o^ y j^j^^ _ 



(19) 



Now for = M$o (i-e. when the magnetic flux is an integer number M of flux 



quanta $o 



2e> 



the condition 



5,T 



holds and we can simultaneously diagonalize 

the operators H (x, y), S, T . By introducing adimensional coordinates on the torus, Eqs. 
(19) can be rewritten as: 



Sip {u) = (a; + 1) , Tip {u) = e^^^'^^^ + t) , 



(20) 



where uj — x + ry, x E [0,1], y E [0,1]. One can look for eigenfunctions 99(0;) = 
giTTMy '^f(^u;) and define magnetic translation operators Sa,Ta acting only on the corre- 
sponding holomorphic part f (cu): 



SJ{u;)=f{u; + a), TJ{i^)= e-^(«^-+2a.) ^ ^ 
In this way Eqs. (20) become: 

Sip (a;) = e'^^y^^S^f (a;) , Tip (a;) = e'^'^y'^TJ (u) . 



(21) 



(22) 



It is straightforward to show that the ground state wave function for M Cooper pairs is 
a coherent state given by: 



iPa (u^i, ...,uJm) = e-^-S-i^?/„ (^1, ...,u;m) , 



(23) 



M 



fa (a;i, ...,ijJM) = Yi 

i<j=l 



Oi (^»i,^) 



Xa (^|t) , 



(24) 
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where ui = X^j^^ (^i is the "center of charge" variable and the non local functions Xa ("^k) 
are the characters of our theory (the TM), whose analytical expressions are reported in 
the Appendix. In fact it can be shown that such characters are eigenfunctions of the 
following generalized magnetic translations operators: 

— Ili=l^a/M^ '^a — Yli=lX^/M ^ (25) 

where 5*^/^ and T^,jj\j are the magnetic translation operators for the single Cooper pair. 
In particular, by recalling the explicit expressions of the conformal blocks (48)- (52) and 
(44)-(47) in the four topological sectors [3], given in the Appendix, and the transformation 
properties of the Theta functions [39], we get (for the step a — 1, see also Eq. (22)): 

^x:|(wc|t) = x+(wc|r), 
P sector, 

'TxI+){wc\t) =ixl+){wc\T), Tx^_^{wc\t) = -ixl_){wc\T), (27) 
A sector, where we defined X(±)('"^c|'?') = X(o)('"^c|'?') ± x'^^^{wc\t), 

^xJ))K|r) = xj^(wc|r), TxJ)(we|r) =x5)(^i^c|r), (28) 
P sector, and 

^X(o)(^ck) = ^X(o)(w^ck), Tx^i)(wc|r) = iX(i)(w^ck), (29) 

for the A— A sector respectively. From such a picture it is evident then how the degeneracy 
of the non perturbative ground states is closely related to the number of primary states 
and that our characters represent highly non local functions: all the topological properties 
of our system are codified in such functions. 

Let us now focus on the general topological structure of the vacua on the torus and on 
their corresponding physical JJA realization in relation to the various boundary conditions 
imposed by closing the square lattice. Let A be the cycle of the torus which surrounds the 
hole of the Corbino disk and B the cycle in the radial direction. On a pure topological 
ground we expect for the torus a doubling of the ground state degeneracy, which can be 
seen at the level of the conformal blocks (characters) of our TM (see Appendix). Indeed 
we get for periodic boundary conditions (i. e. an even number of plaquettes) along the 
74-cycle an untwisted sector, P — P and P — A, characterized by periodic (P) and twisted 
(A) boundary conditions along the S-cycle respectively, described by the four conformal 
blocks (48)- (51); likewise we get for twisted boundary conditions (i. e. an odd number 
of plaquettes) along the A-cycle a twisted sector, A — P and A — A, characterized by 
periodic and twisted boundary conditions along the i?-cycle respectively, described by the 
four conformal blocks (44)- (47). The four different topological sectors just obtained can 



for the P — 



for the P - 



for the A — 
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(a) (b) 

Figure 2: The two vacua in the P — A sector. A straight domain wall is generated along 
the A— cycle. 

be put in correspondence with four minimal configurations for the closed fully frustrated 
JJA as follows. 

Regarding the untwisted sector, in Fig. 1 we show the two vacua corresponding to the 
P — P sector, characterized by a JJA with = 4 plaquettes along the circular direction 
and N = 2 plaquettes along the radial direction, 4 and 2 being the minimum number 
of plaquettes respectively needed in order to fulfill the even-even boundary conditions. 
Then in Fig. 2 the two vacua in the P — A sector, implemented by a JJA with even-odd 
boundary conditions along the two cycles respectively, are sketched. 

Regarding the twisted sector, we get for the A — P and the A — A sectors the config- 
urations shown in Figs. 3 and 4 respectively, characterized by the required odd-even and 
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Figure 4: The two vacua in the A — A sector. A straight domain wall is generated along 
both the y4— cycle and the -B— cycle. 

odd-odd boundary conditions. 

In the next Section the relevant two states of the qubit will be found and will lead 
us to propose the configuration shown in Fig. 3 as the minimal one required to build a 
reliable qubit device protected from decoherence. 
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5 Stability and transformation properties of the ground 
states: the physical qubit 



In this Section we will study the stability and transformation properties of the four ground 
states of the JJA arranged in the Corbino disk geometry under an adiabatic elementary 
flux change (±|f ) through the central hole. Because of the finite energy gap to fractionally 
charged excitation states (in complete analogy with the presence of a gap separating the 
ground state from higher energy states in the Laughlin Hall fluid [40]), such an adiabatic 
transformation is believed to leave the system in a ground state which can be different 
from the original one, due to the occurrence of the ground state degeneracy. This analysis 
will lead us to the identification of the two states of a possible protected qubit and to the 
definition of a tool for its manipulation. 

We use the results provided by our TM model in Section 4 and in the Appendix in 
order to analyze such properties by standard conformal techniques. As we showed in the 
previous Section, in the torus topology the characters of the theory are in one to one 
correspondence with the ground states and any one of them describes the component 
of the "center of charge" for the corresponding ground state wave function. On a pure 
topological ground the torus shows a doubling of the ground state degeneracy, which can 
be seen at the level of the conformal blocks (characters) of our TM. Indeed, as outlined 
in Section 4, we get a rich structure for the vacua corresponding to different boundary 
conditions at the ends of the square lattice. That gives rise to an untwisted sector, P — P 
and P — A, described by the four conformal blocks (48)-(51), and a twisted sector. A — P 
and A — A, described by the four conformal blocks (44)-(47) respectively. Now we are 
going to extract from such vacua the two states of the "protected" qubit. 

Let us notice that the ground state wave functions of the twisted and untwisted sectors 
of the TM are characterized by different monodromy properties along the A-cycle. In 
particular the characters of the untwisted sector are single-valued functions along the bi- 
cycle while the characters of the twisted sector pick up a common (—1) phase factor along 
the A-cycle. Such phase factors can be interpreted as Bohm-Aharonov phases generated 
while a Cooper pair is taken along the 74-cycle. The above observation evidences a strong 
difference between the two inequivalent topological sectors, the untwisted and the twisted 
one respectively, on the torus. Indeed in the twisted sector the ground state wave functions 
show a non trivial behavior implying the trapping of a half flux quantum (| (^)) in the 
hole of the Corbino disk. Instead in the untwisted sector, due to the single-valued ground 
state wave functions, only integer numbers of flux quanta can be attached to the hole. 
It is worth pointing out the central role played by the isospin (or neutral) component of 
the TM in producing the discussed non trivial monodromy properties. To this end let us 
recall that the TM is a c = 2 CFT, composed by a c = 1 charged and a c = 1 isospin CFT 
components, as it is well evidenced by the characters decompositions given in Section 
4. Furthermore the transport of a Cooper pair along the A-cycle can be implemented 
by a simultaneous and identical translation Awc — Au^^ = 2 of the charged and the 
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isospin variables. The charged characters have trivial monodromy with respect to this 
transformation, being: 



Ki{w, + 2\t) = Ki{w,\t), / = 0,..,3, (30) 

while the isospin contribution is the one responsible for the non trivial monodromy of the 
complete ground state wave functions: 

Xo,i(2|r) = Xo,i(0|r) , Xi(2|r) = (-l)xi(0|r) (31) 

and the same is true for the characters Xi3- Observe that the change in sign in the 
last relation of Eq. (31) shows the presence in the spectrum of excitations carrying 

fractionalized charge quanta. More precisely the presence in the isospin component of 
one twist-field (with conformal dimension A = 1/16) characterizes all the conformal 
blocks of the twisted sector and accounts for the trapping of a half flux quantum in the 
hole of the Corbino disk. 

We are now ready to study the stability and transformation properties of the ground 
state wave functions when a magnetic flux change takes place through the central hole of 
the closed JJA. The above analysis shows that at the level of the wave functions it has the 
effect to change the monodromy along the A-cycle due to the corresponding change in the 
Bohm-Aharonov phase. Such a modification can be implemented on the center of charge 
component of the wave function, i.e. the characters, with a well defined transformation. 
In the case of the c/iar^'ec? component this analysis has been brought out already in [41] for 
the quantum Hall effect. Let us adapt here the results for the charged component of our 
TM. On a pure physical ground the fact that we are considering a magnetic flux change, 
which is on one side integer in the flux quantum (one flux quantum change ±|f ) and on 
the other side adiabatic suggests both that the monodromy properties do not change and 
that the system remains in a degenerate ground state. Such a physical picture is in fact 
confirmed for the charged component of our TM; indeed, the flux change is implemented 
on the charged characters by the transformation 

ri%i^,(we|r)=e(^)'^-^+'^i^, (we + ^|r) =i^i+i(we|r), / = 0,..,3. (32) 

Let us notice that the charged component wave functions realize a flip process (Z — > Z + 1) 
under one magnetic flux quantum change. 

However the analysis for the complete TM, with charged and isospin components, is 
more involved, due to the non trivial interplay between charged and isospin components 
summarized in the so-called m-ality parity rule, which characterizes the gluing condition 
for the charged and isospin excitations. The main point being the compatibility between 
such parity rule and the transformation of the complete characters of the TM under the 
insertion of a magnetic flux quantum through the hole of the closed JJA, which reads as: 

V/Kkc|r) = e2-("^-+«(-"+-^))/(^^„ + aT\w, + ar|r) , (33) 

a=l/2 
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where f{wn — 0\wc\t) stays for any character of the TM. More in detail the full list of 
character transformations is as follows. 

In the untwisted sector, we have that the two ground state wave functions of the P — A 
sector decouple, being 

Ti/2X^o)(0|we|r) = 0, Ti/2X^,)(0|we|r) = 0. (34) 
Concerning the P — P sector, we have: 

T^/2xU0\w,\r) = (35) 

and 

Ti/2xJ(0|w;e|r) = X+(0|w;e|r) ( T,/2x:^{0\w,\t) = xj(0kc|r) ). (36) 

Such transformations show the instability of the P — P sector under the insertion of a flux 
quantum through the hole of the closed JJA. More precisely the state Xq (0|'^c|t) decouples 
while the state Xsi^l'^clT) gets excited to the state with a kink-antikink configuration 
X^(0|wc|t) (sec Eq. (52) and comments afterwards). 

Furthermore in the twisted sector, we have that the two ground state wave functions 
of the A — A sector decouple, being 

ri/2X^o)(0kc|T) = , Ti/2x'[,){0\wc\t) = 0. (37) 

Concerning the A — P sector, we have that the two ground state wave functions transform 
as: 

Ti/2xJ))(0kc|r) = X("i)(0|wc|r) , Ti/2X^^){0\wc\t) = x5)(0|wc|r) . (38) 

Concluding, the A — P sector only is stable under the insertion of a magnetic flux quantum 
through the central hole, with the two ground states flipping one into the other under an 
adiabatic flux change of ±|^. That allows us to make the following identiflcations: 

|0) - xJ)(Okc|T), |1) - xfi)(Okc|T), (39) 

in terms of the ground states center of charge wave functions (characters). Then |0) 
and |1) are the two ground states of the closed JJA characterized by an odd number of 
plaquettes along the A-cycle and an even number of plaquettes along the S-cycle. If we 
choose a closed path along the A-cycle, we can deflne a topological invariant as the size 
invariant sum over all its plaquettes X^pXp; which results equal to —1 and +1 for |0) and 
|1) respectively. 

The results just obtained lead us to propose the JJA closed in a Corbino disk geometry, 
characterized by an odd number of plaquettes along circular direction and an even number 
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of plaquettes along radial direction, as our protected qubit. In fact the two ground states 
|0) and |1) work as the two logical states of the qubit and the required one qubit operations: 

|0)^|1), |1)^|0), 

are simply implemented by insertion of a flux quantum (±|f) through the central hole. 
That provides a tool for the control and the manipulation of our device, which is indeed 
a "flux" qubit [27]. 

In Fig. 3 we showed the minimal configuration for such a device, that is a closed fully 
frustrated JJA with = 3 plaquettes along the circular direction and N = 2 plaquettes 
along the radial direction, 3 and 2 being the minimum odd and even number of plaquettes 
respectively needed in order to fulfill all the above requests. 

Given the two degenerate ground states (39) and the operation mode defined in (38) 
it should be possible to prepare the qubit in a definite state and to realize all elementary 
one-qubit operations. We showed in particular how to implement the NOT operation, 
the Hadamard gate will be implemented in a future publication. All that is accomplished 
by an adiabatic flux change of through the central hole of the Corbino disk. 



6 Conclusions and outlooks 

The ground states and the topological properties of fully frustrated Josephson Junction 
arrays (JJA) arranged in a Corbino disk geometry for a variety of boundary conditions 
have been investigated in detail, employing a twisted CFT approach [2, 3]. In this way 
we built up a low-energy effective field theory which is a topological one, due to the 
well known relation with Chern-Simons gauge theories in (2 + l)-dimensions [19]. The 
vacua obtained have been classified in four different topological sectors, according to 
the different boundary conditions imposed at the ends of the square lattice, and put in 
correspondence with four minimal configurations for the closed fully frustrated JJA. Then 
a careful analysis of the stability and transformation properties of the conformal blocks 
corresponding to such vacua allowed us to propose a solid state qubit, protected from 
decoherence, whose operation mode is based only on adiabatic magnetic flux changes 
through the central hole of the Corbino disk. 

Josephson junction arrays have been fabricated within the trilayer Nb/Al — AlOx/Nb 
or Pb/Sn/Pb technology as well as all aluminum technology and experimentally inves- 
tigated [42], but in such a case the application of an external transverse magnetic field 
is needed in order to fulfill the requirement of full frustration and that could be another 
source of decoherence. It is now possible to avoid such a problem by realizing arrays with 
a built-in frustration. In fact high-Tc Josephson junction arrays have been recently pro- 
posed [43], which support degenerate spontaneous current states in zero magnetic field 
due to the presence of plaquettes containing an odd number of 7r-junctions [44]. Such 
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unconventional junctions can be realized because of the d-wave symmetry of high-Tp su- 
perconductors [45], which produces a 7r-shift in the phase of the wave function on one 
side of the junction. Furthermore vr-junctions can be obtained also with super conduct ing- 
ferromagnctic-supcrconducting (SFS) [46] and superconducting-insulator-ferromagnctic- 
superconducting (SIPS) [47] structures. One can also achieve the same effect using a 
barrier which effectively flips the spin of a tunneling electron, i. e. when the barrier is 
made of a ferromagnetic insulator [48], of a carbon nanotube [49] or of a quantum dot 
[50] created by gating a semiconducting nanowire. In this way it is possible to avoid the 
external frustration bias but, in any case, external magnetic fields are needed for control 
and read-out operations. So in principle an experimental setup for the realization of our 
protected qubit can be easily conceived, whose basis is a fully frustrated JJA arranged in 
a Corbino disk geometry with an odd number of plaquettes along the inner hole and an 
even number of them along the radial direction. 
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Appendix: the TM on the torus 



In the following we summarize the whole primary field content of our theory, the TM, on 
the torus topology [3]. 

On the torus [3] the TM primary fields are described in terms of the conformal blocks of 
the Z2-invariant c = | sub-theory and of the non invariant c = | Ising model, so refiecting 
the decomposition on the plane outlined in Section 3. The characters Xo{0\t), Xi(0|r), 
XJ_(0|r) express the primary fields content of the Ising model [51] with Neveu-Schwartz 
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(Z2 twisted) boundary conditions [3], while 



X(of/'(0|«;e|r) = xo(0|r)/^oK|r) + Xi(0|r)X2K|T), 
Xl;f\0\w,\r) = x±iO\r){K,{w,\T) + Ks{w,\r)), 
X(2f^'(0|«;e|T) - xi{0\r)Ko(wc\r)+xo(0\r)K2(wc\r) 



(40) 
(41) 
(42) 



represent those of the Z2-invariant c = | CFT. They are given in terms of a "charged" 
Ka{wc\T) contribution: 



K2i+i{w\r) 



e 



2l+i 
4 



(2u;|4r), with / = 0, 1 and i = 0, 1 , 



(43) 
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and a "isospin" one X/3(0|r), (the conformal blocks of the Ising Model), where Wc — 

- — In Zc is the torus variable of the "charged" component while the corresponding argu- 
2m 

ment of the isospin block is = everywhere. 



If we now turn to the whole c = 2 theory, the characters of the twisted sector are given 



by: 



xfo)(OK|T) = x^(0|t)(xo + Xi)(0|t)(Xo + X2)K 
xfi)(OK|T) = x^^m (xo + X^) {0\t) {K, + K,) {w, 

for the A — P sector and by: 

Xlo)i^\wc\T) = x^^{0\t){xo-Xi){0\t){Ko-K2){w, 
X^i)(OK|t) = xx(0|t)(xo-X|)(0|t)(Xi + X3)K 
for the A — A one. Furthermore the characters of the untwisted sector 
r) = (xoXo - XlXi) iO\T)Ko{wc\T) + (xoXl - XlXo) (0 



r), 

r), 

r), 
r), 

are [3]: 

t)K2 {wc\t), 



(44) 
(45) 

(46) 
(47) 



(48) 



(49) 



X(o)Ww. 

X(i)(0|wc|r) = (xoX^ - XiXo) {0\r)Ko{wc\T) + (^XoXo - XiXi) {Q\t)K2 (Wc|r), 

for the P — A sector while for the P — P sector we have: 

X^(Okc|r) = \{xo- X|) (0|t) (xo - X|) (0|t)(Xo - K2){wM) , (50) 
1 / 



X${Q\w,\t) = 2 + ^0 + ^0 + ^2)(wc|r), (51) 



and 



X+(0|we|T) = X4(0k)X4(0|^) (^1 + ^3) (wc|t). 



(52) 



Let us point out that the above factorization expresses the parity selection rule (m-ality), 
which gives a gluing condition for the "charged" and "isospin" excitations. 

Furthermore, in the P — P sector, unlike for the other sectors, modular invariance 
constraint requires the presence of three different characters. The isospin operator content 
of the character x^(0|wc|t) clearly evidences its peculiarity with respect to the other states 
of the periodic (even ladder) case. Indeed it is characterized by two twist fields (A = 1/16) 
in the isospin components. The occurrence of the double twist in the state described by 
X^(0|'u;c|t) is simply the reason why such a state is a periodic state. Indeed, being an 
isospin twist field the representation in the continuum limit of a kink (or equivalently a 
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half flux quantum trapping) , the double twist corresponds to a double half flux quantum 
trapping, i.e. one flux quantum, typical of the periodic configuration. Indeed the state 
described by X7 (0|'^c|t) embeds in the continuum limit a kink-antikink excitation, i.e. it 
represents an excited state in the P — P sector. In this way, as it happens for all the other 
sectors, the P — P sector is left with just two degenerate ground states (Xa (0|'"^c|''") and 
X^(0|wc|t)) and, as expected on a pure topological base, the ground state degeneracy in 
the torus topology is the double of that of the disk. 

Finally the partition function of the TM model on the torus has the following factorized 
form (see [3]): 

Z{w,\t) = Z''''{wc\t)Zj{t), (53) 

in terms of the Moore-Read partition function (c — 3/2) and of the Ising partition 
function Zj (c = 1/2). They have the following expression: 



7MR 



K|r) = |xr(Okc|r)f + |xr(Oke|r)f+ |x2'''^(0|^e|r)f , (54) 



Zjir) = \Xo{r)f + Mr) + Xi(r) . (55) 
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